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Abstract 

The eigenvalue density of a quantum-mechanical system exhibits oscillations, determined by the 
closed orbits of the corresponding classical system; this relationship is simple and strong for waves in 
billiards or on manifolds, but becomes slightly muddy for a Schrodinger equation with a potential, 
where the orbits depend on the energy. In special cases the simplicity has been restored by rescaling 
the size of the orbit, and perhaps a coupling constant at the same time. We point out that the 
goal can be achieved for any system of this class simply by rescaling the overall coupling constant. 
In each of these situations we inspect critically the relation between the oscillation frequency and 
the period of the orbit; in many cases it is observed that a characteristic length of the orbit is a 
better indicator. When these matters are properly understood, the periodic-orbit theory for generic 
quantum systems recovers the clarity and simplicity that it always had for the wave equation in 
a cavity. Finally, we comment on the alleged "paradox" that semiclassical periodic-orbit theory is 
more effective in calculating low energy levels than high ones. 
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I. INTRODUCTION 



The basic principle relating classical periodic orbits to quantum spectra has been suc- 
cinctly stated by Delos and Du [p]]: 

The average density of states as a function of energy [is] equal to a smooth 
monotonic function, related to the volume occupied by the energy-shell in phase- 
space, plus a sum of sinusoidal oscillations. The wavelength and amplitude of 
each oscillation are respectively correlated with the period and the stability of a 
periodic orbit of the system. For given energy resolution AE, only those periodic 
orbits are significant for which the period is less than 2irh/AE. 

An example of a precise theorem along these lines can be found in ||. (Other expository 
works include || | g, |.) 

This principle, after a prehistory associated with the mathematicians Poisson and Selberg, 
sprang into being in its general form approximately 30 years ago, in three research programs 
dealing with three distinct contexts: 

1. Balian and Bloch considered a bounded region in space. The differential operator 
(corresponding to a quantum Hamiltonian H) is the Laplacian with standard (Dirichlet 
or Neumann) boundary conditions. The relevant classical paths are straight lines with 
specular reflection at the boundary. 

2. Gutzwiller || [J studied a Schrodinger equation in infinite space, with a potential 
that is responsible for confining particles to compact regions and creating a discrete 
energy spectrum. The periodic orbits involved are those of the full classical dynamics 
governed by that potential function. (Balian and Bloch also studied this situation in 
a later paper |T0| .) 



3. The paper inspired a sequence of works by mathematicians p\ [12j , [L3| set in a 
compact Riemannian manifold without boundary. Among the very general operators 
studied in the case most pertinent for us is the Laplace-Beltrami operator for 
the manifold plus a potential function that is treated as a perturbation. That is, 
the manifold itself is responsible for confining the orbits and making the spectrum 
discrete, and the periodic orbits that appear in the theorem are the closed geodesies 
of the Riemannian metric, independent of the potential. 

Although ultimately the most important thing about these three bodies of work is their 
similarity, there are some characteristic differences among them. 
First, Gutzwiller works with the Schrodinger equation, 

i^ = H*l> (h = i), (l.i) 

the others (except Colin de Verdiere |llj) with the wave equation, 

(i.2) 

or its first-order pseudodifferential form 



,f = v^. (1.3) 
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(Indeed, the authors of [13| go so far as to say that "No statement like [our theorem] holds 
if we replace the first order operator [V~H] by a higher order operator.") 

Second, Gutzwiller and the many physicists who have followed him usually speak of 
periodicities in the energy spectrum, determined (reciprocally) by the times (periods) of the 
classical orbits. But the conclusions of the other works are formulated in terms of periodicity 
in wave frequency, to (essentially the square root of energy, E) , determined by the lengths of 
the orbits. A casual reader might think that the distinction between energy periodicity and 
frequency periodicity arises from the respective concentration on the Schrodinger equation or 
the wave equation; in other words, that the square root relating E to uj is merely the square 
root relating the quantum Hamiltonian to the pseudodifferential operator y/~H appearing 
in the wave equation (|1.3|) . That would be wrong. Indeed, in hindsight, it is obviously 
wrong: The conclusions of the theory are statements about the spectrum of the same elliptic 
second-order differential operator H. These facts cannot depend upon which time-dependent 
equation, ( |1 . 1|) or (|1.3|) , was used as a technical tool in discovering them. 

A third difference is that the orbits studied in j7| and [11], [12|, |13| are energy-independent 
objects, but the closed orbits in a potential vary with the energy. The latter fact means that 
the parametrization of spectral oscillations by orbits can be meaningful only over rather 
short energy intervals, in general. (An exception is systems with homogeneous potentials, 
where the orbits at different energies are related merely by dilation in space. See Sec. [VI] and 
|5|, text surrounding (1.40)].) Friedrich and Wintgen ||, [14| pointed out that this problem 
can be at least partially avoided by broadening the point of view to allow variation of a 
coupling constant with energy. The primary aim of the present paper is to generalize their 
observation, showing that it has several variants, one that is applicable to any system, and 
others that may be more appropriate for systems with certain scaling symmetries. In each 
case the natural variable of the spectral oscillations is a certain power of the energy. It is 
argued that the identification of the frequencies of these oscillations with the periods of the 
orbits is rather forced; some other quantity, often a length, is more pertinent. 

These matters occupy most of the remaining sections of this paper. Sec. [TV] is a digression 
to establish a (known) relation among action, energy, and period when the coupling constant 
is not scaled; in the context of scaling, this relation no longer holds, but neither is it necessary. 
In Sec. |T] we elucidate the relationship between the wave (|1.2[) - (p..3|) and Schrodinger ( |1 . 1|) 
approaches (and their correlates in classical phase space); this point is understood by workers 
in the field but is seldom spelled out. (To keep the treatments of the two approaches as 
parallel as possible, we discuss the wave equation in the terminology of relativistic quantum 
theory; this should not be allowed to obscure the applicability of periodic-orbit analysis in 
other, more classical, contexts, such as optics and acoustics, where quantization is just a 
metaphor.) Finally, Sec. [IX| addresses the frequently remarked-upon counterintuitive fact 
that the "semiclassical" periodic-orbit method is more effective in reproducing low-lying 
eigenvalues than eigenvalues in the regime of large quantum numbers. 

Bibliographical remarks: The scaling analysis of the diamagnetic Kepler problem has a 
long history, which has been reviewed by Hasegawa et al. Jl5| as well as Friedrich and Wint- 
gen [Q]. Until the (more recent) paper of Creagh and Littlejohn |IB| this work concentrated 
on the separate eigenspaces of axial angular momentum (as is natural from the quantum- 
mechanical point of view), which somewhat complicated and obscured the connection with 
the classical orbits. The paper |J was the inspiration for the present work, after which the 
author became aware of intervening publications by Friedrich [O, [TS| and Main et al. [HS 



which come close to making the main points of the present paper. The main difference is 
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that here we carry out the analysis in terms of simple "active" transformations — replace- 
ments of orbits by different, geometrically similar orbits — whereas in the earlier work the 
arguments seem to involve "passive" redefinitions of physical units. (Moreover, those coor- 
dinate transformations were noncanonical, forcing the introduction of a varying "effective 
Planck constant" — a concept that the present author finds more confusing than helpful.) 
Our reference list is not exhaustive; many other relevant papers can be traced from the ones 
cited here. 



II. WAVE AND SCHRODINGER DYNAMICS 

A Hamiltonian system is defined by a function H (x, p) of x G Q, a (^-dimensional region 
or manifold, and p G R d . Its classical equations of motion are 



Energy is conserved: Each trajectory of the system (ZTT) remains on a phase-space subman- 
ifold • 

H(x,p) = E (2.2) 

of dimension 2c? — 1. 

Suppose that H (x, p) is quadratic in p and nonnegative, and consider a new Hamiltonian 
function 

h^,p) = H( X ,p) 1 / 2 . (2.3) 
Calling the new time parameter r, calculate the resulting equations of motion: 

^ = Vp h = ±#- 1/2 ^ , ^ = \H- 1 ^ . (2.4) 
dr dt dr 1 dt 

It follows that the trajectories of the new system in phase space are exactly the same as 
those of the old one, but parametrized differently. At a fixed energy the two time scales are 
simply related by 

- = 2VE. (2.5) 

The relation between p and the velocity is changed by the same factor. 

Consider the very special case H = p 2 /2m. Then p is a constant vector with magnitude 
p = y/2mE, and the velocity in system ( |2.1| ) satisfies 

v = — = — , v = — , (2.6 
at m m 



while the velocity in system fl2.4Q satisfies 

dx 1 p 2 1 

u = — — = — = — , u — — = constant. (2.7) 
dr 2VE m 2m 

The trajectories of H dynamics in space-time are, of course, those of a nonrelativistic particle: 
straight lines of arbitrary slope depending on energy. The trajectories of h dynamics in space- 
time are independent of energy, as appropriate for a relativistic massless particle. (Clearly, 
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in this context m does not have the physical significance of a mass; rather, it parametrizes 
the wave speed as c = (2m) -1 / 2 .) Passing to a quantum system by the prescription 

p^_, Vx , E ^ % ^ ( 2 - 8 ) 



converts the relativistic energy relation fo(x, p) = E to a special case of ( |1.3|) , which is 
equivalent to the wave equation 



d 2 i\) 1 



V 2 ^ (2.9) 



<9t 2 2m 

supplemented by a condition of positive frequency. The nonrelativistic energy relation ( |2.2|) 
of course becomes the free Schrodinger equation, a case of ( |1.1| ). 

The energy independence of relativistic trajectories leads to a very direct relationship 
between spectral oscillations and the singularities of the fundamental solution of the wave 
equation ( |1.3| ). This is the meaning of the remark of Duistermaat and Guillemin quoted in 
Sec. I. An up-to-date exposition of that methodology is given in . 

Now turn to a Hamiltonian including a potential energy function, 

ff(x,p) = ^ + V(x), (2.10) 

where (for the moment) ^(x) is nonnegative. The momentum is no longer constant, and 
the velocities satisfy 

v 2 = — [E-V(x)] (2.11) 
(the usual nonrelativistic kinetic-energy relation) and 

P 2 

u — 



2m(p 2 + 2raV(x)) 




E 



(2.12) 



The trajectories in space-time are now curves, whose slopes at each point are greater than 
they would have been for a free particle. In particular, in the relativistic case the trajectories 
stay always inside the local light cones, as appropriate for a relativistic particle with mass 
(induced here by an x-dependent interaction). The main point is that the traces of these 
orbits on configuration space are the same as in the nonrelativistic case. 

In passing, note that if V were allowed to be negative, the classical relativistic trajectories 
would be "tachyonic" according to (2.12). Nevertheless, the wave equation (T72) still obeys 



hyperbolic causality. This situation was studied by Schroer and others in the context of 
quantum field theory |2"I |. We shall not pursue this case further, because it interferes with 



the definition of the square roots in ( |2.jj| ) and (|1.5|). 



III. SCALING A BILLIARD 



Henceforth we adopt the nonrelativistic, or Schrodinger-equation, point of view and re- 
store the conventional constant %. We consider Hamiltonians of the form fl2.10p (without 
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restriction on the sign of V). The quantization prescription 
second-order elliptic partial differential operator 



8|) leads unambiguously to a 



(3.1) 



(In more general cases some variant of Weyl quantization must be chosen to resolve operator- 
ordering ambiguities [El 



The complete specification of the problems (|1.1|) - (|1.3|) may also 



require boundary conditions, which we take to be of the standard Dirichlet, Neumann, or 
Robin type. 

It is assumed that the spectrum of the operator H is discrete, at least below some 
threshold. The essence of the Gutzwiller trace theory 0, g |, |, |, |, g |, |, 
that the density of states, p(E), contains a term 



is 



Py(E) = a 7 sin 



Sty 



+ r/ 7 



(3.2) 



for each classical periodic orbit 7 of energy E; here S^(E) is the action § p-dx of the orbit 7, 
and the details of the amplitude a 1 {E) and the phase 7/ 7 do not concern us now. (Also, this 
simplified formulation does not do justice to the complications that arise when the orbits 
are not isolated and unstable. We touch on such matters in Sec. VIII.) 

Let us examine what ( |3.2j ) asserts for a system of the type studied by Balian and Bloch [|7| . 
The classical paths are straight lines in a region Q C R d with specular reflections at the 
boundary of Q. (The distinction among different boundary conditions shows up in the 
phases r? 7 .) In particular, these paths are the same for all E. The speed is fixed at v — 

2E/m (see ( |2.6| )). If the length of such a closed polygonal path, 7, is L 7 , then its period 



is 



T,(E) 



m 

2E °" 



(3.3) 



The momentum is parallel to the path with fixed magnitude p = \j2mE, so the action is 



S 7 (£T) = j p • rfx = V2mE L 7 = 2£T 7 



(3.4) 



Thus 



becomes 



p-y(E) = a 7 sin 



a 7 sin 



2T 7 £ 



+ ^7 



V2mE 



(3.5) 
(3.6) 



It is crucial to note that T 7 is itself a function of E, but L 7 is constant. Thus p 7 in ( |3.6| ) is 
exactly and globally (to the extent that a 7 and r^ 7 can be regarded as constant) a sinusoidal 
function of cu = \J~E~ with period 

P,„ = \l— (3.7) 




2m 



or P w = 2tt/L 7 if 

equation. This is the result obtained by Balian and Bloch M working from the relativistic, 



1, the most natural normalization when dealing with the wave 
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or wave-equation, point of view: The eigenvalue density oscillates with frequency L 7 as a 
function of uo. 

On the other hand, we have 



~dE 



m 

2E^ 



T — ^2- 
7 IE' 



Approximating S^(E) locally by a linear function M, one gets from ( |3.2|) and (|3 



p-y{E) ~ a 7 sin 



(3.8) 



(3.9) 



Thus p 7 oscillates locally as a function of with approximately defined frequency T^{Eq)/% 
(not twice that, as might appear from (|3.5|) ), or period 



2nh 



(3.10) 



This is the sort of statement about spectral oscillations that is most common in the physics 
literature surrounding the Gutzwiller trace formula. 

For completeness we record the time-dependent action 



where 



FLy(T) = S~(E) ETy = ET-y 

mL 2 



mL 2 
2T 



E = \mv 2 



2^2 



has been used. It satisfies 



dRy 

It 



-E, 



(3.11) 

(3.12) 
(3.13) 



which is dual to (3~B) when ( 3.11 ) is interpreted as a Legendre transformation. It can also 
be obtained as 



IL= & Cdt = ET, 



(3.14) 



where £ is the Lagrangian, here equal to the kinetic energy. 



IV. INTERLUDE: THE ACTION PERIOD RELATION 



In converting (|3.2|) to (|3.9| ), only the relation 



dS^y 

~dE 



(4.1) 



was used. This identity is frequently cited as "a well known theorem of classical mechanics", 
but a proof is hard to find in the literature. Indeed, at first it is not even clear what 
( flip means for a generic Hamiltonian ( |2.10 ), since 5 7 was defined for a single closed orbit 7 
existing at a particular energy, say Eq . The discussion in Sec. |T| shows that for a billiard the 
orbit indeed persists unchanged (in configuration space) as E varies; but when the dynamics 
involves a potential function the classical path defined by an initial point and direction on 
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7 will generally cease to be a periodic orbit when E deviates from Eq . Generically 7 will 
smoothly evolve with E into a family of nearby periodic orbits ^{E) (in general, disjoint 
from 'y(Eo)), but even this picture breaks down at certain singular points (e.g., a maximum 
of the potential in dimension 1, where two families of orbits merge into one). The main 
point of the remainder of this paper is that such issues can be evaded by adopting a broader 
point of view. Here, however, we provide a simple derivation of ( |4.1| ) in circumstances where 
it makes sense. 

Assume the existence of a family ^y(E) of classical closed orbits of ( |2.1| ) depending 
smoothly on the energy, E (as E varies over some interval, possibly small). Assume also 
that either (1) the period T 7 depends smoothly and monotonically on E, so that 7 can be 
alternatively parametrized by T, or (2) T 7 is independent of E. In Case (1) we shall prove 

dR.y „ , . , 

-w = ' < 4 - 2 > 

which is equivalent to (|4. 1|) under the Legendre transformation 

S = R + ET (S = S(E), R = R(T)). (4.3) 

For a family of classical trajectories from a fixed initial space-time point (x', 0) to a variable 
final point (x, T) (and satisfying condition (1)), ( |4.2| ) is the usual (e.g., [p3~| ) Hamilton- 
Jacobi equation for R as a function of (x, T) (with E = if(x, Vi?)); but for periodic orbits 
a different argument is necessary. 
By definition, 

R= I p ■ dx- ET = f Cdt, (4.4) 

J-y JO 

C = fmx(t) 2 - V(x(t)). (4.5) 
In the standard way, consider a variation 

5R= / —6xdt + C\ T 6T, (4.6) 
Jo ox 

where 

5C = mx • <5x — W • <5x = — (mx • <5x) — (mx + W) ■ <5x. (4-7) 

dt 



The last term in (|4.7| ) vanishes by the equation of motion, and integrating the other term 
yields 

mx(T) ■ 5x(T) - mx(0) ■ 5x(0) = -mx(T) 2 5T, (4.8) 



because 
so that 



x(T)=x(0), x(T) = x(0), (x + 5x)(T + 5T) = (x + 5x)(0), (4.9) 



5x(T) = (x + 5x)(T) -x(T) 

= (x + 5x)(T + 5T) - (x + <5x)(T)oT - x(0) 

= 5x(0) -x(T)5T + 0((5T) 2 ). (4.10) 
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mx(T) 2 + §mx(T) 2 - V(x(T)) 



Putting all the ingredients into (|4.6| ) yields 

5R 

~5T ~ ~ 

= -imx(T) 2 - V(x(T)) 

as desired. 

In Case (1), (|4.2|) and ( [4.3| ) imply ( [4.1| ) by the standard Legendre argument. In Case (2), 
the Hamilton- Jacobi equation ( |4.2| ) does not exist, but in that case the foregoing calculation 
shows that 5R = 0, and hence ( [4.1[ ) follows from ( |4.3[ ) even more directly. In either case, 



then leads as before |TJ to ( |3.9j ) and ( |3.10| ), under the additional assumption that T, 



7 

is sufficiently slowly varying for the concept of a local frequency of oscillation of p 7 to make 
sense. 

An interesting example (and the principal instance of Case (2)) is the one-dimensional 
harmonic oscillator, 

V(x) = \mw 2 x 2 . (4.11) 
The general solution of its equation of motion is 

x(t) = xq costot + — sinwt. (4-12) 
m 

T = -, E = ^- + \mw 2 xi ) (4.13) 

S = ET='^, R = 0. (4.14) 
w 

dS T= S 1 (4.15) 



One calculates 



Thus 



dE E 
in agreement with (|4.1|) (but in contrast to (|3.8| )). 

V. SCALING A COUPLING CONSTANT 

Consider a Hamiltonian 

ff(x,p) = |l + AV(x) (5.1) 
with a generic potential V and a coupling constant A. The equations of motion are 

* -»W M , (5.2) 
at m at 

and any solution curve lies on an energy surface 

^- + \V{x) = E. (5.3) 

Let (x (t), Po(t)) be a solution of (|5.2|) and ( |5.3|) with A = A and E = E Q . For any 
positive real number a, consider 

x(t) = x (at), p(t) = apo(at). (5.4) 
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A short calculation shows that (x(t),p(£)) satisfies (|5.2|) and ( |5.3|) with 

A = a 2 An , E 



a 2 E 



(5.5) 

In particular, the path x(t) in configuration space is independent of the parameter a 2 = 
E/Eq. In this sense the same closed orbit exists for all values of the energy. When V = 
this is the familiar billiard orbit reviewed in Sec. |TTT| . When the potential is not zero, 
however, one must pay the price of varying A — i.e., considering different physical systems 
at different energies — to gain the convenience of fixed orbits. (However, in electromagnetic 
wave problems with dielectrics the scaling ( |5.5| ) does relate solutions of the same system p4}|.) 
For the orbit ( |5.4| ) it is clear that 



a 



(5.6) 



and one calculates 



S 



j> p • g?x 



ft 



ft 



Po(ftt) 



To 



Pol^J 



Thus (|3.2|) becomes (by use of (|5.5|)) 

P ,(E) : 



a 7 sm 



a 7 sm 



T 




[ftx (ftt)] dt 
x (t) dr = a S . 



(5.7) 



■ UJ 



7? 7 



That is, just as for a billiard, we have a globally sinusoidal function of uj = 
frequency Sq/Tiujq , or period 

2irhuj 



S 



Thus the analogue of the L 1 of a billiard is (cf. (|3.7| )) 

FfF 2n So 



2m R, 



2m\E \ 



(5.8) 
E\ with 

(5.9) 
(5.10) 



For the harmonic oscillator orbit ( 4.12|) , this characteristic length is ir times the maximum 
value oix(t) (not 4 times, which would be the actual length of the orbit). It should be noted 
that (|5.9| ) or (|5.10| ) depends only on the orbit 7, not on the reference scale arbitrarily chosen 
to correspond to ft = 1. 

One can also look at the local frequency with respect to E, following ( p.8| )-( pTT0[ ). From 
fl5.7| ) and ( |5.5| ) one has 

^ = A. (5.11) 
dE 2E K J 

But this quantity is no longer equal to T, because of the variation of the coupling constant 

with E. For the harmonic oscillator it equals |T, by ( |4.15[) . More generally, since the 

momentum is parallel to the velocity for a system of type (|5.l|) , one can compute 



S = j> p • kdt = l£[E ~ AV(x(t))] dt, 



(5.12) 
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(5.13) 



and hence 

S m A 

— = T 

2£ £ 

Note also that S/2E scales as 1/a and therefore cannot be regarded as a characteristic time 
of the orbit 7 (in contrast to the length (]5.10 )). 

Finally, because nearby nonperiodic orbits also obey the scaling law, the amplitude a 7 
in (|5.8|) (determined by the monodromy of the Poincare map) is independent of E, whereas 
the same cannot be said for the a 7 in ( |3.9| ). The phase shift r? 7 also is constant for Dirichlet 
and Neumann boundary conditions (not Robin). 

Main et al. [|19| rescale the kinetic term of the Hamiltonian instead of the potential term 
and the total energy. Their procedure is equivalent to that of this section together with a 
redefinition of the time scale. In that approach the energy of a given orbit is always fixed, 
and the quantum parameter (see Sec. |V111| ) is y/m instead of uj. 



VI. SCALING A HOMOGENEOUS POTENTIAL 



When the potential in ( |5.1|) is homogeneous of degree u, 

V(J3hl) = /TV(x) (for (3 > 0), 



(6.1) 



an alternative scaling procedure exists that avoids changing the coupling constant. Assume 
that the configuration space Q, if not all of R , is invariant under dilations (e.g., a half-space 
or a cone) with boundary conditions that are either Dirichlet or pure Neumann (not Robin) . 
(If Q does not satisfy this condition, then it must be dilated along with the orbits, and again 



the system has been replaced by a one-parameter family.) In place of ( |5.4| ) consider 



x(t) = a 2 Xo(u t), p(t) = a u p (a u - 2 t). 



(6.2) 



This is a solution of (|5.2j ) and (|5.3| ) with 

A = A , 



E = a 2u E . 



(6.3) 



The orbits for different values of a are not the same, but they are geometrically similar; one 
can say that they are the same if the unit of length is rescaled by a 2 . 
We now have (parallel to fl5.6|)-( 



a 



T = a 2 ' u T Q 



Py(E) = a 7 sin 



E , 



S = a u+2 S 



„ f +2 

S n / E\ — 



h \E, 



+ ^7 



(6.4) 
(6.5) 
(6.6) 



Thus the variable with respect to which the spectral oscillations take place is a peculiar power 
of E, and the analogue of the characteristic length ( |5.10| ) is a rather inconvenient function 
of the fiducial action and energy. Note that the exponent of \E\ in (|6.6| ) equals 1 for the 
harmonic oscillator; that it acquires the value | (familiar for billiards) only as v —>■ +00; 
that it is negative for — 2 < v < 0; and that the cases v — and v = — 2 are singular. 
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The analogue of Q5.11D is 



dS v + 2 S 



dE 2v E 
( |5.12| ) remains true, and the analogue of (|5.13|) is 



v + 2 S v + 2 



T-±J Q T V(x(t))dt 



2v E v 

On the other hand, the theorem of Sec. [TV] applies to the orbit family (|6.2j), and so 



(6.7) 



(6.8) 



dS_ = T 
dE 



Equating (|6.8|) and ( |6.9|) yields the identity 



t 2ET 
\V(x{t))dt = — 
o v + 2 



(6.9) 



(6.10) 



for a periodic orbit of a homogeneous potential of degree u, which is a rewriting of the virial 
theorem . 



The Coulomb potential, 



Wx) = -- (r 2 = x 2 + y 2 + z 2 ) } 
r 



has 



v = -1, 



v + 2 1 



(6.11) 



(6.12) 



2v 2 

Furthermore, in the regime of discrete spectrum and periodic orbits of the atom, E is 
negative. Thus ( |6.7| ) becomes 

(6.13) 



dE 2\E\ ' 

which displays a superficial coincidence with the corresponding relation for a billiard, (KBi 



VII. SCALING OF MIXED TYPE 



A mixture of the strategies of Sees. |V| and |V| is appropriate in some situations. Suppose 
that the potential is a sum of two terms, each of which is homogeneous: 



7x(/3x) = fVjlx), V 2 (/3x) = /T 2 ^(x). 
The equations of motion and the energy equation are 

-77 = - , -£ = — AiVVi(x) - A 2 VK 2 (x) 



i^ + A 1 V 1 (x) + A 2 V 2 (x) =£?. 
2m 



(7.1) 

(7.2) 
(7.3) 
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(This is not a perturbative situation: Neither potential is assumed to be small compared 
to the other.) Following the usual pattern, suppose that we have a solution with coupling 
constants Aio and A20 and energy E . Perform the scaling (|6.2|) appropriate to V\ : 



x(t) = a 2 xo(a^" 2 t), p(t) = a Ul p (a Ul - 2 t). 
Then Q solves Q and Q with 

Ai = Ai , A 2 = a 2 ^-^A 20 , 



E = a 2ui E . 



(7.4) 



(7.5) 



Equations (|6.4|)-(|6T6|) continue to apply, with v = v\. 

The point of this transformation is that a fixed orbit can be associated with a one- 
parameter family of classical situations, in which one coupling constant is held fixed but 
the other varies with energy. The prototype is the Hamiltonian studied by Friedrich and 
Wintgen §0, 

(x, p) = f- - - + \mw\x 1 + y 2 ). (7.6) 
2m r 

It describes, in a rotating coordinate system, a hydrogen atom in a constant magnetic field 
B = (2mc/e)tu along the z axis. This is "a real physical system that can be and has 
been studied in the laboratory". The magnetic field is a continuous variable that is under 
the experimenter's control; the charge of the proton is not! Therefore, to get a family of 
experimentally realizable systems one applies (|7.1| ), (|7.4| ), ( |7.5|) with v\ = — 1, Vi = 2. Thus 



A, 



a 



A20 j 
becomes 



a 7 sin 



and the bottom-line equation 

Pl {E) 

In the notation of Friedrich and Wintgen, 

r = 



So 



E 



E 



a 



~ 2 Eq 



'1/2 



W 

w 



a 



E 
Eo 



T 2/3 



(7.7) 



(7.8) 



(7.9) 



and the contribution p 7 to the spectral density is a sinusoidal function of r -1 / 3 . (We write 
T and w for the 7 and uj of || to avoid notational collisions.) For the beautiful fruits of this 
approach as applied to ( [7.6| ) we refer to [b|, |T^] . 



From S = a Ul+2 S n and E 



a: 



2vi 



Eo we still get 



dS 
dE 



m + 2 S 
2~v x E 



(7.10) 



but this is not equal to an orbital period in these incompletely homogeneous systems. 

More generally, let us repeat the analysis of (|7.1|) - (|7.5|) , but with an arbitrary potential 
in the role of V2 ■ Instead of a scaling of A 2 as in ( [7.5|) , one obtains a one-parameter family 
of potentials 

V 2[a] (x) =^y 2 (a- 2 x). (7.11) 

The classical orbits persist as E and x rescale and V 2 deforms in this way, while the ho- 
mogeneous potential V\ and its coupling constant remain unchanged. When V2 is a sum 
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of several homogeneous terms with different exponents Uj, a separate scaling law for each 
coupling constant, 

Xj = a 2( - Ul ~ Uj) X j0 , (7.12) 

results from Q7.11J) . This last situation has been discussed thoroughly by Friedrich |17| , [18| , 
with special attention to the easily experimentally realizable case where V\ is a Coulomb 
potential and the Xj are the strengths of applied constant magnetic and electric fields. 



VIII. IMPLICATIONS AND SUMMARY 

A central feature of the relation between classical and quantum mechanics is that quan- 
tization introduces into each problem a new fundamental scale, set by the quantum of 
action, h. That is, one-parameter families of situations that are equivalent classically be- 
come distinct quantum-mechanically. Taking the semiclassical limit refers to motion along 
one of these families in a certain direction, the opposite direction leading to "deep quantum" 
behavior. Moving from one family to another, on the other hand, corresponds to various 
purely classical distinctions, such as integrable versus chaotic, or different values of angular 
momentum. 

When the potential is homogeneous (Sec. |V|) or absent (Sec. |Tj), a "situation" simply 
means a point in phase space. In the more general contexts of Sees. [V] and |V11| the space of 



situations should be enlarged by one dimension, representing the coupling constant A or A2 or 
the parameter a in ( |7.11| ). In either case the points are grouped into equivalence classes, the 



classical trajectories. Finally, the trajectories (including, in particular, the closed orbits) fall 



into families, ( |5.4| ) or ( |6.2| ), related by a geometrical similarity and therefore mathematically 
equivalent. The classical trajectories of a family are exactly the same at all energies, except 
for a trivial rescaling. But the quantum states at various energies are quite different. 

In the study of spectra the most convenient choice of parameter along each family is 
the energy. What determines whether an energy is large or small? To compare it with 
h one must construct a quantity with dimensions of action. The dimensions of J2m\E\ 



are action divided by length, so a suitable measure is the product of J2m\E\ with some 
length characteristic of the entire system. The statement u h is small" is meaningless unless 
translated into such a criterion. In the end one can always choose units in which h = 1 and 
2m = 1, and then the only independent physical dimension is length. (Once the quantum of 
action and the characteristic length of the system are fixed, rescaling m amounts to changing 
the unit of time. Since m has been totally inert throughout our considerations, it could have 
been eliminated at the very beginning, but so that the classical-mechanical equations would 
look familiar, we did not do so.) The criterion now is whether the length li^ -1 / 2 = uj~ 1 
is large or small relative to the length scale of the system. This restatement makes sense 
in classical wave theories, such as optics or acoustics, as originally studied by Balian and 
Bloch 0, m . 



The length scale can be set by the geometry of the boundary, if there is one. Otherwise 
it must be a characteristic of the potential. At first sight this is a simple matter: Quantities 
such as V 2 V/V (evaluated at a minimum of the potential, say) are independent of the 
coupling constant and characterize the spatial scale of the potential, as distinct from its 
strength. However, if V is homogeneous, (|6.1|) shows that the distinction between the spatial 



scale and the coupling constant is a mirage; and Sec. VII suggests that this ambiguity can 



be imported into more general potentials as well. For a homogeneous potential there is an 
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A 



FIG. 1: Solid curves: Energy levels of the harmonic oscillator, E n = (n + i^jfiw, as functions of 
coupling constant, A oc w 2 . No quantum states exist in the shaded area. Dashed lines: Loci of 
typical classical orbits fixed in space (Sec. |V]). Dotted lines: Loci of typical homogeneously scaled 
orbits (Sec. VI). The direction of increasing E along either set of lines leads into the semiclassical 
region of large quantum numbers. 



alternative way to set the scale of the system, exploiting ( pTT]) to "transmute" the coupling 
constant into a length, xq = f3~ l . For example, if d = 1 and the potential energy is Xx u , 
then A has dimensions [energy] [length] ~ u and can be written 



A = Xq 

2m 



(u+2) 



so that 



\x h 



K 1 



2mxr 2 



(8.1) 



(8.2) 



; " \x , 

It is instructive to look at everyone's favorite homogeneous potentials, the harmonic 
oscillator ( |4.11| ) and the hydrogen atom (|6.11|) , in the light of the foregoing remarks. If we 
treat these systems according to Sec. [V|, we are taking a fixed "classical situation" to be 
a fixed ratio E/X, with A oc w 2 or e 2 respectively. We expect to encounter semiclassical 
behavior as \E\ — > oo along one of the diagonal lines in Figs. |T] and fj, and indeed this leads 
into the region of large quantum numbers in either case. There is no unique way to associate 
a natural length with either of these potentials, but lengths independent of A can be built 
out of derivatives of V at an arbitrarily chosen point. If we treat the systems according to 



Sec. [VI], we keep A constant and therefore move along one of the vertical lines in the figures. 
From 



we have 



x Q 



Xq 



1 2 1 

— , E n x = n + - 
w I 



(8.4) 
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FIG. 2: Solid curves: Energy levels of the hydrogen atom, E n oc — e 4 /ra 2 , as functions of coupling 
constant, A oc e 2 . No quantum states exist in the shaded area. Dashed lines: Loci of typical 
classical orbits fixed in space (Sec. |V|). Dotted lines: Loci of typical homogeneously scaled orbits 
(Sec. VI). The semiclassical region of large quantum numbers is reached by moving along the 
dashed lines in the direction of increasing \E\ {E — > — oo), or along the dotted lines in the direction 
of decreasing \E\ (E — > 0). 



for the energy levels of the oscillator, and 

1 



x£ = e 4 , E n xo oc - ^ (8.5) 



rr 



for those of the atom. The semiclassical regime of large quantum number n is reached by 
going vertically upward in either of the figures. For the oscillator this is again the limit of 
large E, albeit in a different direction in the E-X plane. But for the Coulomb potential 
this limit corresponds to E — > 0, not \E\ — > oo. (In view of ( |6.13|) this is still the direction 



of increasing S.) This (known) result is not surprising in view of the negative exponent in 
(|6.6|), ( |6.13| ); but it indicates that the identification of "semiclassical" with "large energy" 
can be ambiguous, even within the theory of the same physical system. 
Let us summarize the main points: 

1. The concept of "spectral oscillations" with respect to energy (associated with the basic 
formula ( |3.2j ), which we took as given) is problematical, because the classical orbit 
structure itself varies with the energy, quantitatively and sometimes qualitatively. 

2. The orbits can be "nailed down" by scaling the coupling constant along with E. From 
this enlarged perspective a mechanical system (with varying A) has one (very clearly 
distinguished) "quantum dimension" along which very well-defined spectral oscillations 
take place. It is necessary to study an entire family of systems (parametrized by A) 
to gain a clear and complete picture of a single system; this is usually not a problem 
theoretically, though it may be very hard experimentally. 
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3. For a homogeneous or partially homogeneous potential there are other options for 
fixing the orbit structure, involving dilations. 

4. The wave frequency u = E 1 ^ 2 is a better variable than E for describing the oscilla- 
tions globally. (In the homogeneous case the exponent \ is replaced by {y + 1)/2za) 
Correspondingly, the best descriptor of the oscillations caused by a particular orbit is 
a characteristic length (5.10) of the orbit, not the orbital period. 



5. Even the local frequency of spectral oscillation with respect to energy, dS/dE, is not 
equal to the orbital period when the coupling constant is varied. (In passing we 
supplied a proof of the equality when A is constant.) 

Points 1 and 2 require further comment in view of the fact that, strictly speaking, the 
Gutzwiller formula (|3.2|) applies only to an orbit that is both isolated and unstable. Similar 



formulas apply when the orbits form continuous families because of symmetries [[Rj, |27 



with completely integrable systems [^8| as the extreme case. Those formulas have the same 
simple oscillatory nature as (|3.2|) , but with different phases and amplitudes; our observations 
apply to them unchanged. (Indeed, the diamagnetic Kepler problem (|7.6j ) is of this type, 
when it has not been reduced to an angular momentum eigenspace and a two-dimensional 
classical configuration space.) Complications develop, however, when symmetric systems are 



slightly perturbed |29], |3(], [31]] or when stable orbits develop bifurcations [|32], |33|, |34]] . In these 
cases the semiclassical formulas involve factors (typically Bessel functions) that interpolate 
between different oscillatory regimes. These phenomena occur even in billiard systems, as 
pointed out in Appendix D] and exemplified in detail in various studies of the ellipse 
f29| , |3T| , |35|j . Therefore, the observations of the present work provide no escape from them. 
The point is, however, that whatever takes the place of ( |3.2j ) must scale in the same simple 
way when the energy and coupling constant are simultaneously scaled as in Sec. V. This 
is so because the semiclassical formulas are built entirely out of classical ingredients, along 
with strategically placed factors of h to give all quantities the correct dimensions. 



Two examples are Sieber's formula for the ellipse [31], eq. (64)] and the same author's 
generic formula for bifurcations [|34], eq. (37)]. In the former case, given a fixed ellipse the 
only relevant parameter is uo = \J~E~ (k, in the notation of [[31]). The formula describes a 
rigid function of uo just like our (|5.8| ), albeit of a more complicated functional form. In the 
other case, in contrast, two things are going on simultaneously. The energy is changing, but 
also the function being evaluated upon the energy is changing. (Most spectacularly, some 
terms completely disappear when the energy passes through the bifurcation point, where 
two "satellite orbits" disappear from the classical structure.) A bifurcation also occurs in 
the billiard problem, but one must vary the eccentricity of the ellipse to see it. By scaling 
the coupling constant along with the energy, one introduces a "normal coordinate system" 
into the manifold of problems under study, so that the quantum parameter uj is decoupled 
from the classical parameters defining the system. In the ellipse problem this was automatic: 
the original energy is already a good quantum parameter, and the eccentricity is the most 



important classical parameter. For a Schrodinger problem ( |5.1| ) the normal coordinates are 
less natural: the quantum parameter is the joint magnitude of E and A; varying it leaves the 
system at a fixed effective "distance" from the bifurcation; and letting A deviate from the 
scale of E amounts to changing the classical mechanics toward or away from the bifurcation. 
Although this approach will not by itself play the major role in deriving more formulas like 
(37)], we submit that it can aid considerably in their interpretation. 
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IX. A FINAL COMMENT 



When periodic-orbit theory was first applied numerically to the spectra of concrete sys- 
tems, early authors expressed surprise that it gave accurate results "even" for the lowest 
eigenvalues. As a semiclassical method, the technique had been expected to be applicable 
primarily in the regime of "large quantum numbers" and hence high energy Later it be- 
came clear that the method was generally practical only in the low-energy regime, because 
the number of periodic orbits increases roughly exponentially with period, and as energy 
increases it becomes necessary to consider increasingly long orbits in order to resolve indi- 
vidual eigenvalues J|, |5[. Long after this fact has been accepted, it is still often regarded as 
"paradoxical". 

We suggest that, in hindsight, this phenomenon is an instance of the familiar principle 
that when a calculation is very stable and involves some kind of smoothing or averaging, the 
inverse calculation is likely to be unstable (highly dependent on the details of the input) and 
hence difficult. (The following remarks obviously lack technical precision and are offered as 
preliminary ideas only. They do not explain how formulas such as (|3.2| ) apparently have 
a wider validity than the stationary-phase arguments that lead to them.) There is a loose 
analogy with the solution of an initial-value problem for the heat equation. For large positive 
time the problem is very easy to solve and very insensitive to the details of the initial data; 
consequently, reconstructing the data from the final solution is hopeless. For very small time 
an approximate solution to the backwards problem, adequate for some practical purpose, 
may be feasible; of course, the forward problem is less trivial in that case. 

In spectral asymptotics we are interested in deducing the spectrum of a differential op- 
erator from its geometry, or vice versa. Here "geometry" is meant in an extended sense, 
including not only the literal geometry of the region Q where the wave functions are de- 
fined, but also the potential function (or other coefficient functions in the operator H). The 
classical periodic orbits of the Hamiltonian i/(x, p) are also aspects of the geometry of the 
system. 

In the "old" spectral asymptotics, associated with the names Weyl, Thomas-Fermi, Mi- 
nakshisundaram, Schwinger-DeWitt, etc., one associates the (high-i?) asymptotic behavior 
of the density of states with the global geometry of the operator: the volume of Q, the inte- 
grated curvature of its boundary, the integral of the Ricci curvature of Q if it is a manifold, 
the integral of the potential V over Q, and so on. Via the asymptotics of the heat kernel, 
the passage from the spectrum to the geometry is rigorously asymptotic. In the inverse 
direction, however, the geometry does not determine a genuine asymptotic expansion of the 
eigenvalue density in powers of E~ l , precisely because of the presence of the oscillatory terms 
that are the subject of the periodic-orbit theory. The formal expansion of the eigenvalue 
density becomes literally asymptotic only when some kind of averaging is performed, such 



as Lorentzian smoothing [26] or Riesz means [36 



A complementary situation has always existed at the low end of the spectrum. If we 
know Q and H exactly, then it is relatively easy to construct the lowest-lying eigenvalues 
and eigenf unctions, by variational methods, numerical methods, etc. Larger eigenvalues 
are harder. Furthermore, merely from a knowledge of the lowest energies one would not 
expect to be able to reconstruct Q or H . (Complete knowledge of a single eigeiafunction is 
a different matter, however.) 

The "new" spectral asymptotics of periodic orbits extends this picture. Decreasing the 
width over which the eigenvalue distribution is averaged, one supplements the power-law 
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asymptotics with the longest-wavelength oscillatory components (|3.2|) , characterized by the 
lengths or periods of the shortest classical orbits. This is still an asymptotic (high-energy) 
matter. It is best thought of as a prediction of classical behavior (well-defined orbits for 
wave packets) in the high-energy, large-action regime. It still involves an averaging of the 
spectrum, albeit on a smaller scale. The stable direction of prediction is from the spectrum 
to the classical, geometrical description. 

The counterpart of this at low energy is that from the classical orbits (and their associ- 
ated amplitudes and phases, cf. (|3.2| )) one can predict the formation of discrete eigenvalues 
(resonant frequencies), with greatest precision near the bottom of the spectrum. In some 
sense this involves an averaging over geometrical information. More precise geometrical 
information in principle allows precise construction of longer-period orbits and hence more 
complete reconstruction of the spectrum (both improved resolution and extension to higher 
energies). The stable direction of prediction is from the geometry to the spectrum. (In the 
other direction, "[Mjissing a level has disastrous consequences for locations, heights, and 
even the very appearance of the non-Newtonian peaks in the Fourier transform." pTfl ) 

It is misleading to think of the regime of low energy (or small quantum numbers) as the 
"deep quantum regime", as if the duality between the spectrum of quantum eigenvalues and 
the spectrum of classical periodic orbits is irrelevant there. A better phrase is "resonant 
regime": There, resonant behavior emerges on the spectral side, out of the geometry and 
classical mechanics, much as, in the opposite "classical regime", classical-mechanical behavior 
emerges on the geometrical side, out of the quantum substrate. 
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